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Abstract. A cycle basis in an undirected graph is a minimal set of 
simple cycles whose symmetric differences include all Eulerian subgraphs 
of the given graph. We define a rooted cycle basis to be a cycle basis 
in which all cycles contain a specified root edge, and we investigate the 
algorithmic problem of constructing rooted cycle bases. We show that a 
given graph has a rooted cycle basis if and only if the root edge belongs 
to its 2-core and the 2-core is 2-vertex-connected, and that constructing 
such a basis can be performed efficiently. We show that in an unweighted 
or positively weighted graph, it is possible to find the minimum weight 
rooted cycle basis in polynomial time. Additionally, we show that it is 
NP-complete to find a fundamental rooted cycle basis (a rooted cycle basis 
in which each cycle is formed by combining paths in a fixed spanning tree 
with a single additional edge) but that the problem can be solved by a 
fixed-parameter-tractable algorithm when parameterized by clique-width. 


1 Introduction 


A cycle basis of an undirected graph is a set of cycles such that all cycles in the 
graph have a unique representation as an algebraic sum of basis cycles. In this 
paper we study algorithms for finding a special type of cycle basis which we call 
a rooted cycle basis , in which all cycles in the basis contain a specified root edge. 


Cycle bases have diverse applications including subway system scheduling 16 
the analysis of distributed algorithms [2], and bioinformatics 1,15 . The specific 


motivation for our rooted variant of the problem comes from mechanical engineer¬ 
ing, where cycle bases have long been used in static analysis of structures such 
as truss bridges 11 and in the kinematics of moving bodies 13 . We recently 


used this method as part of a system for constructing the configuration space 
of moving linkages 19 , systems that include automobile suspensions, fold-out 


sofa-beds, and legs for walking robots. 

In this configuration space construction problem, systems of rigid two-dim¬ 
ensional links are connected at joints where one link can rotate around a point 
of another with one degree of freedom. A system of links and joints is called a 
kinematic chain ; fixing the position of one ground link results in a system called 
a mechanism or inversion , and distinguishing a second input link (connected to 
the ground by a joint and to which force is applied to control the rest of the 


system) results in a system called a linkage 19 . The structure of a linkage can be 











Fig. 1 . A linkage and its linkage graph (a subdivision of K 3,3) 


expressed combinatorially by a linkage graph, an undirected graph with a vertex 
for each link and an edge for each joint, including a distinguished ground-input 
edge. The requirement that the combined motion of the linkage have one degree 
of freedom can be e xpre ssed combinatorially by the property that the linkage 
graph is (|,2 )-tight 14 : every A:-vertex induced subgraph must have at most 


— 2 edges, and the whole graph must have exactly |n — 2 edges, where n is 
the number of vertices in the graph and links in the linkage. Links may cross 
each other in the plane, resulting in a non-planar linkage graph (Figure 1). 

Given a linkage with its linkage graph, each input-to-ground path has an 
associated equation representing the requirement that the joints along the path 
have angles consistent with the fixed ground position at both ends of the path. Our 
system for constructing the configuration space of a linkage chooses a complete 
and non-redundant subset of these path equations and uses Dixon determinants to 
solve this system of equations 19 . Each path determining one of these equations 
can be turned into a cycle by adding the input-ground edge, and a set of equations 
chosen in this way is complete and non-redundant if and only if the corresponding 
set of cycles forms a cycle basis of the linkage graph. However, all of these cycles 
contain the input-ground edge, so the system of equations that we seek comes 
from a rooted cycle basis. Additionally, because the equation solver forms the 
computational bottleneck of our system, we would like the system of equations 
that we construct to be as simple as possible, corresponding to the problem of 
finding a minimum rooted cycle basis. 


New results. We provide the first algorithmic study of the problem of constructing 
rooted cycle bases. We have the following new results: 

— As a warm-up to our main result, we show that an arbitrary graph G with 
designated root edge e has a rooted cycle basis in which all cycles contain e 
if and only if the 2-core of G is 2-vertex-connected and contains e. When a 
rooted cycle basis exists, it can be constructed in time 0{mn). This is tight: 
there exist graphs for which every rooted cycle basis has total size 0{mn). 

— Our main result is that, in an unweighted or positively weighted graph with 
a designated root edge, we can find the minimum weight rooted cycle basis 
by a randomized algorithm with nearly-optimal 0{mn + n 2 logn) expected 
time or by a polylogarithmically slower deterministic algorithm. This basis 
is always weakly fundamental: its cycles can be ordered so that each cycle 
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contains an edge that is not in any earlier cycle. Our algorithm uses a greedy 
method for finding each cycle, with a tie-breaking rule that avoids greedy 
choices that do not lead to a valid cycle basis. 

— In Appendix [C] we show that it is NP-complete to determine whether a graph 
G with root edge e has a fundamental rooted cycle basis, a rooted cycle basis 
determined from a spanning tree T by choosing all cycles formed by an edge 
not in T and a path in T. It remains NP-complete even when G is planar. 
Our proof is based on the observation that, in planar graphs, fundamental 
rooted cycle bases are dual to a form of Hamiltonian cycle. Additionally, we 
use Courcelle’s theorem to show that finding a fundamental rooted cycle 
basis is fixed-parameter-tractable in the clique-width of the input. 


In comparison, for arbitrary cycle bases, every graph has a fundamental cycle 
basis, which may be constructed using any spanning tree algorithm. Finding 
unrestricted minimum weight cycle bases takes polynomial time [2,10,12j |18] . 
However, finding an unrestricted minimum weight weakly fundamental cycle basis 
is NP-hard 21 , and cannot be solved by the same greedy strategy that we use 


for rooted cycle bases, of choosing the shortest cycle that includes a new edge. 


2 Preliminaries 


By F 2 we mean the field with two elements 0 and 1 under mod-2 arithmetic. If U 
is an arbitrary finite set, the subsets of U form a vector space over F 2 with 
the empty set as origin and the symmetric difference of sets as addition. 

We define a rooted graph to be an undirected graph G = (V, E ) with a 
designated root edge e. A cycle is a connected 2-regular subgraph; a cycle is 
rooted if it contains e, and Hamiltonian if it contains every vertex of G. The edge 
space of G is the vector space Ff\ The cycle space of G is the subspace of the 
edge space generated by edge sets of cycles; its elements are subgraphs of G with 
even degree at every vertex 24 . A cycle basis of G is a set of cycles that forms a 
basis of the cycle space 12 . A cycle basis is rooted if all its cycles are rooted. 

A spanning tree of an undirected graph G is a subgraph that includes all 
vertices of G, and is connected with no cycles. Any edge / that does not belong 
to a spanning tree T gives rise to a fundamental cycle for T consisting of / plus 
the unique path in T connecting the endpoints of /. The fundamental cycles for 
T form a cycle basis; a basis formed in this way is called fundamental. 

A matroid 25 may be defined as a family of subsets of a finite set, called the 
independent sets of the matroid, with two properties: 


— Every subset of an independent set is independent. 

— If 1 1 and I 2 are independent sets and |/i| < | 1, then there exists an element 

x belonging to I 2 \ h such that I\ U {a;} is independent. 

The linearly independent subsets of a finite family of vectors in any vector space 
form a linear matroid. In a matroid, a basis is an independent set all of whose 
supersets are dependent; for linear matroids, this notion coincides with the 
standard definition of a basis of a vector space. 
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If the elements of a matroid are given real-valued weights, then the basis with 
minimum total weight can be constructed by a greedy algorithm, generalizing 
Kruskal’s algorithm for minimum spanning trees: initialize a set B to be the 
empty set, and consider the elements in sorted order by their weights, adding 
each element to B if the result would remain independent. In particular, if the 
edges of an undirected graph G are given weights, the weight of a cycle may 
be defined as the sum of the weights of its edges, and the weight of a cycle 
basis may be defined as the sum of the weights of its cycles. Then the minimum 
weight cycle basis may be found by considering all of the cycles of the graph in 
sorted order by weight, adding each one to the basis if the result would remain 
independent. This algorithm may be sped up by considering only a special set of 
polynomially-many candidate cycles, leading to polynomial-time construction of 
the minimum weight cycle basis in any graph j2j[l0,12,18 . 

A simple path in a graph G is a connected subgraph with two degree-one 
vertices (its endpoints) and with all remaining vertices (its interior vertices) 
having degree exactly two. An open ear decomposition of G is a collection of 
simple paths Pi for i = 0,1,2,... (called ears) with the following properties: 


— The first ear P 0 is a single edge. 

— The two endpoints of each ear Pi with * > 0 appear in earlier-numbered ears. 

— No interior vertex of an ear appears in any earlier ear. 


A graph has an open ear decomposition if and only if it is 2-vertex-connected 
(no vertex deletion can disconnect the remaining graph) [26 . This decomposition 
can be constructed in linear time, with any edge as the first ear 17 20 23 . The 
number of ears equals one plus the dimension of the cycle space. 

A vertex of G belongs to at least one cycle of G if and only if it belongs to the 
2-core of G, the subgraph formed by removing isolated vertices and degree-one 
vertices until all remaining vertices have degree > 2. Therefore, the cycle bases 
of G are the same as the cycle bases of its 2-core. 


3 Existence and construction of rooted cycle bases 

The following lemma is a special case of Menger’s theorem, but we give a proof 
as we use the proof construction in our algorithms. 

Lemma 1. Let e be an edge of a 2-vertex-connected graph G. Then for every 
two distinct vertices u and v of G there exist two vertex-disjoint paths (possibly 
of length zero) from u and v respectively to the two endpoints of e. 

Proof. Let Pq = e,Pi,... ,Pk be an open ear decomposition of G. We apply 
induction on fc, with the following cases: 

— As a base case, if k = 0, we have two length-zero paths, one for each endpoint. 

— If k > 0 and neither u nor v is an interior vertex of Pk , the result follows by 
induction on the union of the ears up to Pk- 1 - 
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— If k > 0 and exactly one of u or v is an interior vertex of Pk, without loss 
of generality (by swapping u and v if necessary) we may assume that u is 
the interior vertex. At least one endpoint of Pis a vertex w distinct from v. 
By induction, v and w can be connected by vertex-disjoint paths to e, using 
only vertices in ears Pq, ..., Pk-i- The result follows by augmenting the path 
from w with the part of path Pk from u to w. 

— If k > 0 and both u and v are interior vertices of Pk, then u and v have two 
disjoint paths within Pk to the endpoints of Pk . By induction, the endpoints 
of Pk can be connected by paths to e, using only vertices in ears Pq,..., Pk- 1 - 
The result follows by concatenating these paths with the paths within Pk- 

Thus, in all cases, the desired two paths exist. □ 

An ear with one edge cannot be part of a path constructed by this proof. So 
for a graph G with n vertices and m edges and a known ear decomposition, we 
can discard the one-edge ears and transform the case analysis of the proof into 
an algorithm that constructs the two desired paths in time 0(n). 

Theorem 1. An undirected graph G rooted at edge e has a cycle basis that is 
rooted at e if and only if e belongs to the 2-core of G and the 2-core is 2-vertex- 
connected. When a rooted cycle basis exists, it can be constructed in time 0{mn) 
and the total length of the cycles in the basis is 0{mn). 

Proof. If G has a rooted cycle basis, its 2-core must be 2-vertex-connected. For, 
suppose that a vertex v is deleted from the 2-core. Every remaining vertex u 
belongs to a basis cycle from which only v can have been deleted, leaving a path 
connecting u to the remaining endpoints of e. In this way any two remaining 
vertices can be connected to each other via e, so the remaining vertices are not 
disconnected. 

In the other direction, suppose that the 2-core of G contains e and is 2- 
connected. Then it has an open ear decomposition Pq = e, Pi ,..., Pk- We may 
form a set of cycles C\ , C' 2 ,..., in which each cycle C. t consists of e, the edges 
in Pj, and two paths through the union of ears P l7 P 2 ,... P,_i connecting the 
endpoints of Pi to the endpoints of e. These cycles are independent because each 
one contains at least one edge in Pj that does not belong to any previous cycle. 
As an independent set of cycles of the correct cardinality to be a basis, they must 
be a basis. 

After computing the ear decomposition, each cycle takes time 0{n ) to con¬ 
struct (by the remarks following lemma 1) and has length 0(n), giving the stated 
time and length bounds. □ 


The length and time bounds of the theorem are tight in the worst case: for a 
graph consisting of two 0(m)-vertex cliques connected by two 0(n)-vertex paths 
(Figure 2 left), every cycle through e and an edge in the farthest clique from 
e has length n, so every rooted cycle basis has total length 0(mn). For linkage 
graphs with m = |n — 2, the time and length bounds become 0{n 2 ), which is 
again tight: every rooted cycle basis of an n -vertex ladder graph (Figure 2 right) 
has total length 0(n 2 ). 
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Fig. 2. Two graphs whose rooted cycle bases all have large total length: two cliques 
connected by two long paths (left), and a ladder graph (right). 


In contrast, unrooted cycle bases may be significantly smaller. Every graph 
with m vertices and n edges has an (unrooted) cycle basis of total length 
0(min(n 2 ,rnlogn)), a bound that is close to tight because of the existence 
of sparse graphs of high girth for which every cycle basis has total length 
f2(n log n) [6}[l2]. 

4 Finding the minimum weight rooted cycle basis 

In this section we show how to find a rooted cycle basis of minimum total length 
in biconnected graphs with positive edge weights, in polynomial time. We use 
a greedy algorithm that chooses one cycle at a time, and prove it correct by 
showing that the sequence of cycles selected by this algorithm correspond to an 
ear decomposition. Our strategy is to show that an optimal basis can be derived 
from an ear decomposition: the cycles of the basis can be sorted from shorter to 
longer cycles in such a way that, in each successive cycle, the edges that do not 
belong to earlier cycles form an ear. Our algorithm performs the following steps: 

1. Initialize what will eventually become a cycle basis to the empty set. 

2. Use Suurballe’s algorithm (reviewed in Appendix [A]) to compute, for each 
edge, the shortest rooted cycle through that edge. 

3. While there exists an edge that is not included in any of the cycles chosen 
so far, select an edge that has not yet been included and whose computed 
shortest-cycle length is as small as possible, and add its cycle to the basis. 

We will prove this algorithm correct under the additional assumption that no 
two paths, and no two cycles, have the same weight as each other. We say that a 
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Fig. 3. An unweighted rooted graph (left) with two shortest rooted cycles that together 
cover the whole graph but do not generate its cycle space (center and right). Our 
algorithm requires that no two cycles have equal weight, to prevent bad sets of cycles 
such as these from being chosen. 
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graph is unambiguously weighted when this is the case. When paths and cycles 
can have equal weights, this algorithm can fail by choosing a set of cycles that 


together cover all edges but do not generate the whole cycle space (Figure 3), 
so we need a consistent tie-breaking rule in this case, which we describe in 
Appendix [B] 


4.1 Greedy cycle sequences 

We define a greedy cycle sequence to be a sequence of cycles that could be 
produced by the algorithm described at the beginning of this section. That is, it 
is a sequence of rooted cycles C\, C 2 , ■ ■ ■ in which 

1. Each cycle includes an edge that is not in any earlier cycle in the sequence, 

and 

2. Subject to constraint (1), each cycle is as short as possible. 

We will prove a sequence of lemmas about greedy cycle sequences, with the goal 
of showing that the set of new edges added by each cycle forms an ear and therefore 
that our greedy algorithm for rooted cycle bases is correct. To do so, it is helpful 
to have a notation for the subgraph of G formed by the vertices and edges in the 
first i cycles in the sequence. We call this subgraph the itli ambit of the cycle basis, 
denoted A,-. 

Lemma 2. Let G be an unambiguously-weighted 
rooted graph, Ai be the ith ambit of a greedy cycle 
sequence C\, C 2 , ■ ■ ■ for G, and x be a point in 
Ai (a vertex or a point interior to an edge). Then 
Ai contains the shortest path in G from x to each 
endpoint of the root edge of G. 

Proof. Let t\ and £2 be the endpoints of the root 
edge. We will show by induction on i that A, t con¬ 
tains the shortest path from x to fq; by symmetry 
it also contains a path to £ 2 - We may assume that x 
does not belong to A,_i, for otherwise the shortest 
path is already contained in Aj_ 1 by induction. 

Then let P be the shortest path in G from x to 
ti; we claim that P must remain within Ci until it 
reaches a vertex of A,_|. For, if P deviated from 
Ci at some vertex y outside of A,_i, let z be the 
first point at which P returns to a vertex of Cp, z 
must exist, because P eventually reaches t\, which 
belongs to Ci. In this case the rooted cycle formed 
from Ci by removing the path in Ci from y to z and replacing it with the part of 
P from y to z would be strictly shorter than Ci (because the part of P from y 
to z is a shortest path and no two paths have equal length) and would contain 



Fig. 4. Notation for |Lcmma 2| 
A shortest path is shown di¬ 
verging from Ci before reach¬ 
ing A;_i, proven impossible by 
the lemma. 
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the vertex y outside Aj_i, contradicting the greedy choice of C,; as the shortest 
rooted cycle not contained in A,_i. 

Therefore Ai contains the portion of P from x to i, and by induction it 
contains as well the rest of the shortest path from x to t\. □ 


4.2 Rungs of the Suurballe ladder 

Let Pi and P 2 be two disjoint paths from s (an arbitrary vertex in the given 
graph G) to t\ and t 2 (the endpoints of the root edge of G), as constructed by 
Suurballe’s algorithm. Recall that this algorithm constructs two different paths 
from s to ti and f 2 , the first of which is a shortest path in G and the second of 
which is the shortest path in a derived graph H. The union of these two paths 
differs from Pi and P 2 by a collection of paths that we call rungs. Each rung is 
traversed in one direction by the shortest path in G and in the opposite direction 
by the shortest path in H. The endpoints of the rungs lie on Pi and P 2 (in 
the same order on both paths) and each of the two shortest paths is formed by 
following one of Pi or P 2 until reaching the endpoint of a rung, then traversing 
that rung and continuing in the same way along P 2 or Pi until the next rung, etc. 


Lemma 3. With s, P\, and P 2 as above, let C be 
the cycle formed by Pi, P 2 , and the root edge, and 
let R be the rung of P\ and P 2 closest to s. Let 
D be the cycle formed by removing the parts of P\ 
and P 2 from s to the endpoints of R, and replacing 
them with R. Then C is longer than D. 

Proof. Let u and v be the endpoints of R , let £ su 
denote the length of the path in Pi from s to u, 
let £ sv denote the length of the path in P 2 from 
s to v, and let £ uv denote the length of the rung. 
Pi follows C from s to u then crosses rung R; by 
construction, it is the shortest path in G from s to 
t\. Thus, £ su + i vv < £ sv , for if not then Pi couldn’t 
be a shortest path. Equivalently, adding £ su — £ uv 
to both sides of the inequality gives 2£ su < £ su + 
£sv ~ £uv But the left hand side of this inequality 
is positive (by the assumption that the input graph 
has positive edge weights) and the right hand side 
is the difference in weights between C and D. □ 



Fig. 5. Notation for |Lemma 3| 
Pi and P 2 are shown in con¬ 
trasting colors; the horizontal 
segments with both colors are 
the rungs. 


Corollary 1. Let Ci be a cycle in the greedy cycle sequence, and let A t _-| be the 
ambit of the previous cycle. Suppose that Ci is constructed by applying Suurballe’s 
algorithm from a starting vertex s, and suppose that the two disjoint paths Pi 
and P 2 comprising Ci have a nonempty sequence of rungs. Then these rungs, and 
all parts of Pi and P 2 from the first rung endpoint to ti and f 2 , belong to 1 . 












Fig. 6. Notation for |Lemma 4| The figure shows a cycle Ci such that Ci \ Ai -1 
more than one path, a configuration that is proven impossible by the lemma. 


forms 


Proof. Otherwise the cycle D in the statement of Lemma .1) or one of the other 
cycles constructed in the same way from one of the other rungs, would be a 
shorter cycle containing at least one edge that is not in i, and would have 
been selected in place of Ci in the greedy cycle sequence. □ 


4.3 From cycle sequences to ear decompositions 

As we now show, a greedy cycle sequence with cycles C\, C 2 , • • • and ambits A \, 
A 2 , ... may be used to derive an ear decomposition, in which the first ear Pq is 
the root edge and each subsequent ear Pi is the subgraph Ai \ Ai -That is, for 
each i , this subgraph is a single path. 


Lemma 4. Let Ci and Ai be a greedy cycle sequence and the corresponding 
sequence of ambits in an unambiguously-weighted graph. Then Ai \ Ai_i forms a 
single connected path in the given graph G. 


Proof. Let ti and t 2 be the endpoints of the root edge. As with any rooted 
cycle not contained in A, : _ j, C,; forms one or more connected paths in Ai \ A t _j. 
separated by vertices or edges of A;_i. Let P be the closest to t\ of these paths 
according to their ordering along Ci, and let / and v be the farthest edge and 
vertex from t\ in the path ordering of P. Then v is an endpoint of / and of path 
P, and belongs to A % -\. By Lemma 2] the shortest path from v to t\ stays within 
Ai_i and therefore does not use edge /. 

Let s be a point on edge /, sufficiently close to v that the shortest path from 
s to 1 1 passes through v. (In other words, subdivide the edge, and place a vertex 
at this point.) If we apply Suurballe’s algorithm starting from the point s, it 
will find two paths Q\ (the shortest path from s to ii) and Q 2 - The symmetric 
difference of these two paths is C, ; again, because C t is the unique shortest rooted 
cycle containing s. Q 1 passes from s through v and then stays within Ai- 1 by 


Lemma 2 Q 2 follows the rest of P, then stays within Ai- 1 until it reaches the 


endpoint of the first rung of Q 1 fl Q 2 (by the choice of P as the first of the paths 
in Ai \ A^ 1 in cycle order). After reaching this rung endpoint, Q 2 continues to 
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stay within Aj_i by Corollary 1 Thus, neither Q\ nor Q 2 can escape Aj_ 1 once 
they enter it, so the cycle C t that they form can only have the single component 
P outside of A i _ 1 . □ 


Corollary 2. Let Ci be the cycles of a greedy cycle sequence for the rooted graph 
G, and let Aj be the corresponding ambits. Let Pq be the one-edge path formed by 
the root edge ofG, and for i > 0 let Pi be the path Ai\Ai_\. Then the sequence of 
paths Pq, Pi,... Pi is an ear decomposition of the subgraph Aj, and the sequence 
of all paths formed in this way is an ear decomposition of G. 

Proof. By |Lemma 4[ each of these graphs is a path; its endpoints belong to earlier 
paths and its edges and interior vertices do not. Thus, this sequence of paths 
satisfies all the requirements of an ear decomposition. □ 


4.4 Greed is Good 


We have nearly completed the proof of correctness of our greedy algorithm for 
constructing minimum weight rooted cycle bases. 


Lemma 5. Let Ci be the cycles of a greedy cycle sequence for an unambiguously- 
weighted rooted graph G, and let Aj be the corresponding ambits. Then the set of 
cycles Ci, C 2 , ...Ci is a minimum weight rooted cycle basis for Aj. 


Proof. We use induction on i. For i = 1, the graph C\ = A\ has only the one 
cycle. For i > 1, this set of cycles is linearly independent because each contains 
at least one edge not found in earlier cycles. The number of cycles is the same as 


the number of ears (after the root edge) in an ear decomposition, by Corollary 2 


so it equals the dimension of the cycle space. As an independent set of the correct 
number of cycles, these cycles must form a cycle basis for Aj. 

Because the cycle space of any graph forms a matroid, the minimum weight 
basis of any subset S of cycles can be found by a greedy algorithm that at each 
step selects the minimum-weight member of S that is independent of previous 
selections. By the induction hypothesis, any cycle that is independent of the 
previous i — 1 selections must use at least one edge outside of Aj_i, and Cj is 
the minimum-weight rooted cycle with this property. Therefore, the cycles form 
a minimum-weight rooted cycle basis. □ 


Theorem 2. The minimum weight rooted cycle basis of a biconnected rooted 
graph G with positive edge weights can be constructed in polynomial time. 

Proof. As outlined at the beginning of this section, we use Suurballe’s algorithm 
to order the edges of G by the lengths of their shortest cycles through the base 
edge. Then, using this order as a guide, we construct a greedy cycle sequence 
by repeatedly choosing an edge / that is not part of the already-chosen cycles 
and using another instance of Suurballe’s algorithm to find the shortest cycle 
through / and the root edge, breaking ties in favor of cycles that use as few new 
edges as possible. By |Lemma 5[ the resulting set of cycles will form a minimum 
weight rooted cycle basis. 
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In a graph with n vertices and to edges, Suurballe’s algorithm can be im¬ 
plemented in time 0(m + nlogn). The first stage of the algorithm may be 
implemented using Dijkstra’s algorithm in this time bound. The second stage 
involves shortest paths in a graph H with negative edge weights, to which Dijk¬ 
stra’s algorithm does not directly apply. However, in this second stage, we may 
re-weight each directed edge in H from u to v with length l, giving it the new 
weight t + d(s, v) — d(s, u), where s is the starting vertex of the first path and 
d is the shortest-path distance between two vertices in the input graph. This 
reweighting does not modify the comparison between any two path lengths, so 
the shortest paths in the reweighted version of H remain unchanged. With these 
weights, the edges whose weights were negative become zero-weight, and all other 
edge weights remain non-negative, so Dijkstra’s algorithm may again be applied. 

A naive implementation of the algorithm applies Suurballe’s algorithm 0(m ) 
times so its total time is 0(m 2 + mn logn). However, this can be improved 
by observing that there are only 0(n) choices for the first path in Suurballe’s 
algorithm, and that for each first path it is possible to handle all starting vertices 
of the second path, simultaneously, by using Dijkstra’s algorithm to perform a 
single-destination shortest path computation. With this improvement the total 
runtime is 0(mn + n 2 logn). The algorithm as described so far applies only to 
unambiguously-weighted graphs but in Appendix [B] we describe how to reduce 
the general problem to this case in polynomial time. □ 
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A Suurballe’s algorithm 


As it is most commonly described, Suurballe’s algorithm 22] (a special case of the 
shortest-augmenting-path method for minimum cost flows) finds two edge-disjoint 
shortest paths between the same two vertices in a weighted graph. However, it 
can readily be adapted to find shortest vertex-disjoint paths from a single source 
vertex to two different vertices, to find the shortest rooted cycle containing a 
given vertex, or to find the shortest rooted cycle containing a given edge. As our 
proof depends on the details of this algorithm, we review them briefly here. 



I 



Fig. 7. Suurballe’s algorithm: a shortest path P\ from s to ti, together with a shortest 
path P 2 from s to 1 2 in a modified graph in which backwards travel along Pi has 
negative cost, combine to form two disjoint paths of minimum total length from s to t\ 
and £ 2 - 


We suppose that we are given a weighted directed graph G with no negatively 
weighted cycles, a designated source vertex s, and two designated destination 
vertices t\ and f 2 - Suurballe’s algorithm begins by finding a shortest path Pi 
from s to t\ (in the graph G\{t 2 }, so that path Pi avoids vertex f 2 )- It then 
constructs a new graph H, as follows: 

— For each vertex v of G , there are two vertices Vi n and u out in H. 

— If vertex v does not belong to path Pi , then H includes a weight-zero directed 
edge from vi n to w ou t • If v does belong to path Pi , then instead H includes a 
weight-zero edge in the other direction, from t> out to V[ n . 

— For each directed edge uv of G that does not participate in path Pi, there is 
an edge with the same weight in H from u out to v- ln . 

— For each directed edge uv of G that does participate in path Pi, there is an 
edge with the negation of its weight in H from w,j n to w ou t ■ 
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There can be no negative cycles in H, for any such cycle would either translate 
to a negative cycle in G or could be used to improve path P 1 . Next, Suurballe’s 
algorithm computes a shortest path in H from s to £ 2 - In terms of the original 
graph G, this corresponds to a path P 2 from s to £2 that travels forwards along 
edges disjoint from Pi, backwards along edges of P\, and cannot touch a vertex 
of Pi without traveling along an edge incident to that vertex. Additionally (with 
the additional assumption that G has no zero-length cycles) the maximal subsets 
of contiguous edges of P\ used by P 2 appear in the same order in both paths, 
for a violation of this ordering would provide a cycle along which one of the two 
paths could be improved. 

The symmetric difference of the two paths Pl and P 2 (the set of edges used 
in one but not both of them) forms the union of the two desired vertex-disjoint 
paths from s to h and t 2 , as shown in |Figurc 7| 

To find the shortest path from ti to t 2 through an edge e rather than a source 
vertex s, we may apply the same algorithm, choosing a shortest path from one 
endpoint of e to ii in G and then either from the other endpoint to t 2 (if the 
first path does not pass through e) or from the same endpoint (if it does pass 
through e). 
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B Making the weights unambiguous 

The basic version of our minimum-weight basis algorithm, described in the main 
text of the paper, requires a graph that is unambiguously weighted: no two paths or 
cycles have the same total weight. We describe here two techniques to simulate an 
unambiguous weighting, for a graph that is not already unambiguously weighted: 
one randomized without any slowdown and one that is deterministic but slower. 

We assume a model of computation in which, however the given edge weights 
are represented, comparisons are exact: any comparison of two sums of weights 
(of paths or cycles in the graph) is computed correctly and deterministically. We 
assume that sums of polynomial numbers of input values can be computed and 
compared exactly, and that our model of computation is capable of representing 
integers with 0(log n) bits of precision (else how could it address the memory cells 
required to store the input) but we do not allow unlimited-precision arithmetic 
beyond these assumptions. 

Our task is to augment the weight comparisons of the cycle basis algorithm 
with extra information so that, whenever two different sets of edges have the same 
weight, we can nevertheless determine an unambiguous perturbed comparison 
result. This perturbation should not change any weight comparisons that were 
already unequal in the given input, it should be performed efficiently, and it 
should not generate any inconsistencies (such as cyclic order relations) that could 
cause our algorithms to produce invalid results. 

B.l Randomized 

Note that, to perform our algorithm correctly, we do not actually need a truly 
unambiguous weighting. It is sufficient to achieve a weighting in which all com¬ 
parisons of weight sums performed by the algorithm are unambiguous. For, if 
we find a perturbed set of weights w for which all algorithmic comparisons are 
unequal, then there also exists a perturbation of the perturbation, w', such that 
w' is truly unambiguous and such that the algorithm makes the same choices on 
w' as it makes on w. But, on w', the algorithm produces a correct output, so it 
must also be correct on w. 

Randomly perturbing the weights of a graph in such a way that our algorithm 
sees only unambiguous weight comparisons, with high probability, is relatively 
straightforward. For each edge e, we choose a random clogn-bit integer r e for 
an appropriate constant c. Then, whenever we compute a sum of edge weights, 
we also compute along with it the sum of these random integers. If the algorithm 
compares two sums of weights and finds them to be equal, it then compares the 
corresponding sums of the ry, and uses the result of that comparison. If the sums 
of the r, are also equal, then we choose arbitrarily which of the two values to 
consider as smaller. 

Any individual comparison of the weights of two distinct edge sets has 
probability at most 2 -clogn = n~ c of finding equal sums of r*. Since the algorithm 
makes at most 0(n 3 ) comparisons in total (in the case of dense graphs), any 
choice of c > 3 will cause the algorithm to succeed with high probability, by 
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the union bound. After the algorithm has run to completion, we verify that the 
result has the property that the newly-covered edge sets of each cycle form an 
ear decomposition with the correct number of ears. If so, it must be a valid 
minimum-weight rooted cycle cover, and if not, a failure must have occurred so 
we choose a new random perturbation and try again. 

We summarize the results of this section: 

Theorem 3. We can find a minimum-weight cycle basis for graphs with n 
vertices and m edges in randomized expected time 0(mn + n 2 logn). 


B.2 Deterministic 

An unambiguous perturbation can always be achieved by choosing an arbitrary 
distinct positive integer index i for each edge, choosing a sufficiently small value 
of <5 (smaller than the minimum difference between the weights of two edge 
sets, and adding 5/ 2* to the weight of an edge with index i. In this way, all 
perturbations are small enough and no two perturbations are equal. It is not 
necessary with this approach to determine an explicitly numeric value for d; 
instead, one can compare the unperturbed weights of edge sets and, if the result 
is equal, compare the perturbation values, as we did in the randomized algorithm. 
However, computing the perturbation values numerically would seem to require 
n-bit numerical precision, not allowed in our computational model. Instead, we 
describe data structural techniques that can determine the same comparison 
values, deterministically, in polylogarithmic time per comparison. 

We define the incremental first difference problem as follows. We have an 
ordered collection of elements (in the rooted cycle basis problem, these are the 
edges of the input graph) and we wish to form sets of elements, starting from the 
empty set, either by adding one element at a time or by specifying the elements 
that belong to a set (in the rooted cycle basis problem, these sets are paths and 
cycles in the graph). Additionally, we wish to be able to compare two sets to 
determine whether they are equal or unequal, and if they are unequal, to find 
the element at which they differ that is earliest in the element ordering. 

To solve this problem, we fix a balanced binary tree T whose leaf nodes are 
the elements and whose leaf ordering is the same as the ordering of the elements. 
For node j of the tree, we let Lj denote the left child of j (another node or a 
leaf element), we let Rj denote the right child of j, and we let Dj denote the set 
of leaf elements descending from j. 

Our representation will consist of a collection of subset objects , each of which 
is associated with one particular node j of T and represents a nonempty subset 
of Dj. If a subset object x is associated with tree node j and represents the set 
X, then x stores pointers for two child subset objects: the left child f x , associated 
with tree node Lj having set X (~l D r. , and the right child r x associated with 
tree node Rj having set X (~l Dr . If either of these two intersections is empty 
then the corresponding child pointers i x or r x may be set to null pointers. In 
this way, one of the given sets in the incremental first difference problem may be 
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represented directly as a subset object associated with the root of T, and by the 
binary tree of subset objects reachable from that root object. 

In order to compare sets efficiently, we require that no two subset objects 
associated with the same node j can have the same associated set. If two of the 
given sets have the same intersection with Dj , then they both share the same 
subset object for that tree node. Thus, we need to be able to find and re-use 
existing subset objects when necessary. To do so, we use a deterministic dictionary 
data structure that we call the subset dictionary. Its keys are the triples (j, l x ,r x ) 
associated with each subset object, and the value associated with such a key is 
the subset object having that triple of values. The subset dictionary could be 
implemented as a binary search tree, with logarithmic cost per operation; however, 
for integer keys with polynomial range, more efficient (in theory) deterministic 
dictionary structures are available with time 0((loglogn) 2 / log log log n) per 
operation 3|. To transform a triple (j,£ x ,r x ) into an integer with polynomial 
range, we assign each subset object a unique index number, we let q be a 
polynomial bound on the number of subset objects to be created, and we calculate 
an integer key as jq 2 + aq + b where a and b are the index numbers of t x and r x . 

An example of the shared subset objects for three sets is illustrated in |Figurc~8l 

Lemma 6. Given two sets with this representation, it is possible to determine 
the smallest element at which the two sets differ in time O(logn). 

Proof. We begin with the two subset objects x and y associated with the root 
node of T that represent the two sets. If these are the same object, then the two 
sets are the same; otherwise, they must either have differing left or right child 
pointers. If the left children differ, we replace x by £ x and y by t y , and otherwise 
we replace x by r x and y by r y , in either case moving down one level in T to 
another pair of differing sets containing the same smallest difference. Eventually 
we will reach a state where either x or y is a null pointer, after which we can 
follow the leftmost non-null path downward from the remaining non-null subset 
object to find the smallest differing element. This process takes constant time 
per level of T so its total time is O(logn). □ 



Fig. 8 . The binary tree T of ordered elements and the subset objects of an incremental 
first difference data structure for the three sets {a, b, d, /}, {6, d, g , h}, and {a, b, /, g}. 
The blue lines show the connections between each tree node or subset object and its 
two children. The set labels on the subset objects are not part of the data structure. 
The subset dictionary is not shown. 
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Lemma 7. Given a set S with this representation, and an element e, we can 
add S U {e} to the representation in time O(logn(loglogn) 2 /logloglog?i). 

Proof. More generally, we show that for any subset object x , associated with a 
tree node j and with associated set X, we can find or construct a subset object 
for X U {e}. To do so, we determine whether e is a left or right descendant of j , 
and let y be the left or right child of x, associated with a set Y. We recursively 
find or construct a subset object y' for the set Y U {e}, then look up the triple of 
j, y', and the other child of x in the subset dictionary. If this lookup succeeds, we 
have found the subset object we desire. If it fails, we create a new subset object. 
In this way, whenever we create a new object we can be sure that its associated 
tree node and set do not coincide with any subset object that has already been 
created. 

This construction process takes a constant number of dictionary operations 
per level of T, or a logarithmic number of dictionary operations overall overall. 
□ 


Using the same bottom-up construction process, looking up each new triple 
in the subset dictionary prior to constructing a new subset object for that triple, 
we can show: 


Lemma 8. Given a set S as a sorted list of its elements, and an incremental 
first difference data structure on a universe of n elements, we can include S in 
the sets represented by the data structure in time O(?r(loglogn) 2 /logloglogn). 

Using this data structure to disambiguate equal weight sums, we have: 


Theorem 4. We can find a minimum-weight cycle basis for graphs with n ver¬ 
tices and m edges deterministically in time 0((mn log n(log log n) 2 / log log log n + 
n 2 log 2 n). 


Proof. We build an incremental first difference data structure on the edge sets 
of the paths and cycles created by the algorithm for unambiguously-weighted 
graphs. This data structure slows down the instances of Dijkstra’s and Suurballe’s 
algorithms within the overall algorithm by the time bound given in|Lemma 7 


as 


these algorithms build up their paths by one edge at a time allowing Lemma 7| 
to be used. Once Suurballe’s algorithm has found two paths whose symmetric 
difference forms a minimum-weight cycle, we can construct a set for that cycle 
using |Lemma 8j 

Within these algorithms, whenever a weight comparison determines that 
two edge sets have the same weight sum, we use |Lemma ~6| to find the edge 
making the dominant contribution to the perturbations for those edge sets, and 
use this information to break the tie appropriately. This modification slows 
each comparison down by a logarithmic factor. In particular, maintaining the 
priority queue in Dijkstra’s algorithm only involves weight comparisons, and not 
construction of any new edge sets, so the 0(n 2 logn) term of our unambiguously- 
weighted algorithm is only slowed down by this factor rather than by the slightly 
larger factor of the 0(rnn) term. □ 
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C Fundamental rooted cycle bases 

A fundamental rooted cycle basis is a rooted cycle basis that is generated by a 
spanning tree, in the sense that all cycles in the basis are formed by a path in 
the spanning tree together with a single non-tree edge. As we show, it is NP-hard 
but fixed-parameter tractable to determine whether such a basis exists. 

C.l Hardness 

We define a plane graph to be a connected graph embedded in the plane with 
no edge crossings; that is, it is a planar graph together with a choice of how 
to embed the graph. The dual graph of a plane graph G has a vertex for each 
bounded or unbounded face of the embedding of G ; each edge in G separates 
two faces (or possibly a single face from itself) and corresponds to a dual edge 
connecting the corresponding two vertices. The dual graph is also plane but not 
necessarily simple: it may have multiple adjacencies (more than one edge between 
the same two vertices) and self-loops (an edge with both endpoints equal to 
each other). It has a planar embedding for which each dual vertex belongs to 
its corresponding face, and each edge of G is crossed once by its corresponding 
dual edge and does not intersect any other dual edge. The dual graph of the dual 
graph of G is isomorphic to G itself. 

We define a rooted Hamiltonian cycle of a graph G rooted at an edge e to be 
a Hamiltonian cycle of G that includes edge e. 

Lemma 9. Let G be a plane graph rooted at e, and let T be a spanning tree of 
G containing e. Then the fundamental cycle basis for T is rooted at e if and only 
if the set of edges dual to the edges in (G\T) U {e} forms a Hamiltonian cycle 
for the dual graph of G. 


Proof. Let V be the set of edges dual to the edges of G\T; then V is necessarily a 
spanning tree of the dual graph of G, and (T, T') form a tree-cotree decomposition 



Fig. 9. A spanning tree generating a rooted fundamental cycle basis of a planar graph, 
and its dual Hamiltonian cycle 
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of G 17,8]. Let P be the path in T' connecting the two faces that are separated 
by e, and let e' be the edge dual to e. 

If P = T ', then T' Uje'} is a Hamiltonian cycle. By the Jordan curve theorem, 
this cycle separates the plane into two parts, an inside and an outside. Let A and 
B be the sets of vertices of G inside and outside the cycle, respectively. Then e 
and every edge of G\T crosses the cycle exactly once, so it has one endpoint 
in A and the other endpoint in B. On the other hand, the edges of T \ {e} do 
not cross the cycle at all, so both of their endpoints lie in the same set. Every 
fundamental cycle induced in T by an edge f £ G\T crosses from A to B at /, 
and must return to A on another edge, but the only other edge that passes from 
one set to another is e. Therefore, every fundamental cycle contains e, and the 
fundamental cycle basis is rooted. 

If P ^ V, then there exists a dual edge /' € T' that is adjacent to a vertex 
v of P but does not itself belong to P. T' U {e'} is not a Hamiltonian cycle, 
because it has degree three or more at v. Again, by the Jordan theorem, the 
cycle P U {e'} separates the vertices into two subsets A and B, with only e and 
the edges dual to P having endpoints in both subsets. Let / be the edge whose 
dual is f. Then the fundamental cycle induced in T by / cannot include e, 
because if it did it would have only one edge with endpoints in both subsets, an 
impossibility. Therefore, the fundamental cycle basis of T includes a cycle that 
does not contain e and the basis is not rooted. □ 

Corollary 3. Let G be a plane graph rooted at e. Then G has a rooted funda¬ 
mental cycle basis if and only if its dual graph, rooted at the dual edge e! of e, 
has a rooted Hamiltonian cycle. 



Fig. 10. If a degree-three vertex in a graph G is replaced by the subgraph on the right, 
then G has a Hamiltonian cycle if and only if the modified graph has a Hamiltonian 
cycle containing e. 


Lemma 10. Let G be a graph containing a degree-three vertex, and replace this 
vertex by the four-vertex subgraph shown in \Figure 10\ Then G has a Hamiltonian 
cycle if and only if the modified graph has a Hamiltonian cycle that contains e. 

Proof. If G has a Hamiltonian cycle, it must use two of the three edges at 
the replaced vertex v. A simple case analysis shows that in each of the three 
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ways of choosing two edges, there exists a corresponding path through the four- 
vertex subgraph that uses the corresponding two edges, covers all vertices, and 
uses e. Therefore, if G has a Hamiltonian cycle, the modified graph also has a 
Hamiltonian cycle that uses e. In the other direction, if the modified graph has a 
Hamiltonian cycle that uses e, then contracting the four-vertex subgraph to a 
single vertex produces a Hamiltonian cycle in G. □ 

Lemma 11. The problem of determining whether a given 3-connected rooted 
planar graph contains a rooted Hamiltonian cycle is NP -complete. 


Proof. We prove the result by means of a reduction from the known N P-complete 
problem of finding a Hamiltonian cycle in a 3-regular 3-vertex-connected planar 
graph 19 . Given such a graph G, modify it according to Lemma 10, giving a 
graph H and an edge e such that we wish to find a Hamiltonian cycle in H 
through e. G has a Hamiltonian cycle if and only if H has a Hamiltonian cycle 
through e, and the transformation preserves planarity and 3-connectivity. □ 


Theorem 5. It is NP -complete to determine whether a given rooted graph has a 
rooted fundamental cycle basis. 


Proof. We prove the result by means of a reduction from finding a Hamiltonian 
cycle through a given edge e of a 3-connected planar graph G (Lemma 111. 
Embed G planarly, let H be its dual graph, and let e! be the edge dual to e. 
Because G is 3-vertex-connected, H is a simple graph. Then by |Corollary 3[ H 
has a fundamental cycle basis rooted at e! if and only if G has a Hamiltonian 
cycle through e. Since finding a rooted Hamiltonian cycle is NP-hard, so is 
testing the existence of a rooted fundamental cycle basis. On the other hand, 
it is straightforward to verify that a given spanning tree T generates a rooted 
fundamental cycle basis, so testing the existence of a rooted fundamental cycle 
basis belongs to NP. Since it is in NPand is NP-hard, it is NP-complete. □ 


C.2 Fixed-parameter tractability 

We now show that it is fixed-parameter tractable (parameterized by clique-width) 
to test for the existence of a fundamental cycle basis. This means that there is an 
algorithm for solving the problem whose running time is a polynomial of the input 
size (independent of the parameter) multiplied by a (non-polynomial) function of 
the parameter. The clique-width parameter that we use is the minimum number 
of colors needed to construct the graph by a sequence of operations that take 
disjoint unions of colored graphs, add edges between all pairs of vertices matching 
a specified color pair, and recolor all vertices of a given color with a different color. 
Graphs of bounded treewidth also have bounded clique-width, but not necessarily 
vice versa. A monadic second-order (MSOi) formula for a graph property is 
an expression that combines variables representing vertices or sets of vertices, 
quantifiers over these variables, terms that test whether a vertex belongs to a set 
or whether two vertices are adjacent, and the standard Boolean connectives over 
these terms. The connection between these two concepts is given by the following 
variation of Courcelle’s theorem: 
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Lemma 12 (Courcelle, Makowsky, and Rotics |5|). Any graph property 
that can be expressed by an MSOi formula has a fixed-parameter tractable algo¬ 
rithm when parameterized by clique-width. 

Thus, to prove fixed-parameter-tractability, we need only find an MSOi 
formula expressing the existence of a fundamental rooted cycle basis. The following 
lemma accomplishes this. 

Lemma 13. Let G be a graph rooted at edge e. Then G has a fundamental rooted 
cycle basis if and only if there exist two vertex sets Si and S 2 such that all of the 
following conditions are met: 

— Si and S 2 partition the vertices of G: no vertex belongs to both sets and every 
vertex belongs to at least one of the two sets. 

— The two endpoints of e belong to different sets of the partition. 

— For each i, set Si induces a connected set. That is, there does not exist 
another set C such that SiC\C and Si\C are both nonempty and such that 
there are no edges connecting Si D C to Si\C. 

— For each i, set Si is acyclic. This can be expressed by the property that Si 
induces a graph with an empty 2-core: for each nonempty subset D of Si, 
there exists a vertex v G D that has at most one neighbor in D. 

Proof. If G has a rooted fundamental cycle basis generated by tree T, then let 
Ti and T 2 be the two subtrees of T formed by deleting edge e, and let Si and S 2 
be the vertex sets of Tj and T 2 . Then these two sets partition the vertices, and 
induce connected subgraphs. The subgraphs they induce are also acyclic, because 
if the subgraph induced by Si included an edge / that was not part of Ti , then 
/’s cycle in the fundamental cycle basis would not be rooted. 

Conversely, suppose that G and e have sets Si and S 2 that meet these 
properties. Then S\ and S 2 induce two disjoint subtrees of G which together 
span all the vertices of G; adding e to these two subtrees gives a spanning tree T. 
Every edge / in G \ T must connect one of these two subtrees to the other; the 
path in T between the endpoints of / must necessarily pass through e, because 
e is the only edge in T that spans the cut from Si to S 2 . Thus, T generates a 
rooted fundamental cycle basis. □ 

Corollary 4. The existence of a fundamental rooted cycle basis can be expressed 
by an MSOi formula. 

Theorem 6. Determining whether a given rooted graph has a rooted fundamental 
cycle basis is fixed-parameter tractable when parameterized by the clique-width of 
the graph. 

Proof. This follows immediately from|Lemma 12|and|Corollary 4| □ 
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